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at  various  times  after  the  vail  has  been  struck,  and  pressure -time  histories 
of  the  gas  at  the  vail.  Representative  values  for  the  peak  vail  pressure  at 
various  distances  from  the  center,  for  the  constant  gamma  case,  are  22,712 
Bars  at  6.0  cms.,  005.6  Bars  at  24.0  cms.,  and  5 • 37  Bars  at  126.4  cms. 
Scaling  to  other  distances  and  charge  veights  is  considered.  A  512-lb 
pentolite  charge,  for  example,  is  shown  to  produce  a  pressure  of  145  psl  at 
the  wall  some  26  ft  away.  Such  a  pressure,  if  applied  to  a  spacecraft, 
would  probably  destroy  it. 


Approved  by: 

J.  7.  Moulton,  Jr. ,  Chief 
Air-Ground  Explosions  Division 


PUBLISHED  NOVEMBER  1962 


EXPLOSIOSS  RESEARCH  DSPARTMDIT 

U.S.  NAVAL  OREHANCE  LABORATORY 
WHITE  OAK,  MARYLAND 


1 


VOI/St  62-19  15  January  1962 

Explosions  in  Vacuum 


19m  proxies  of  explosions  In  rarefied  atmospheres  can  b«  attacked  by  numerl- 
cal  methods  using  high-speed  electronic  computers.  This  report  deecribea  a 
preliminary  calculation  in  vhich  a  apharlcal  charge  detonatea  in  vacuum, 
eauaing  the  resultant  detonation  gaaea  to  expand  and  strike  a  rigid,  spheri¬ 
cal  target.  The  work  was  done  under  HOL  Task  HUHE-3-E-000/ 212-1/ F008-10- 
004  (PA  0l4) ,  ana  is  a  partial  solution  of  KEY  PROBLEM  Ho.  1$,  Chapter  7* 
IAT0KD  Report  3906,  (Key  Problems  In  Explosives  Research  and  Development) . 


V.  D.  OOUMUI 
Captain,  081 

Commander 


a, 


C.  J.  AROHSOH 
By  direction 


Tk' 


NOLTR  62-19 


CONTENTS 

Page 

INTRODUCTION .  1 

INITIAL  CONDITIONS .  1 

EQUATION  OF  STATE .  2 

DESCRIPTION  OF  THE  NUMERICAL  COMPUTATION .  4 

RESULTS  .  5 

SCALING .  6 

CONCLUDING  REMARKS .  7 

REFERENCES  .  24 

ILLUSTRATIONS 

Figure  Title 

1  General  Form  of  7  Versus  p  Curre  .  3 

2  Wall  Pressure  Versus  TIm  (Wall  at  8.0  Cm.) .  8 

3  Wall  Pressure  Versus  Tins  (Wall  at  16.0  Cm.) .  9 

4  Wall  Pressure  Versus  Tins  (Wall  at  32.8  Cos.) . 10 

5  Wall  Pressure  Versus  TIm  (Wall  at  64.4  Cos.) .  11 

6  Wall  Pressure  Versus  Time  (Wall  at  126.4  Cm.) .  12 

7  Peak  Pressure  Versus  Wall  Distance . 13 

8  Pressure  (Megabars)  Versus  Distance,  Wall  at 

8.0  Cm,  Ganna  -  2.682  . 14 

9  Energy  Density  (Energy/ Volume  x  4*R2,  in  Units  of  Negater  - 
Cm2).  Wall  at  8.0  Cm,  Gamma  ■  2.682  Area  Under  Curve 

Equals  Energy  in  Units  of  Megater  -  Cm3 . 15 

10  Energy  Density  (Energy/ Volume  x  4*R2,  In  Units  of  Megater  - 
Cm2).  Wall  at  8.0  Cm,  Gasan  ■  2.682  Area  Under  Curve 

Equals  Energy  in  Units  of  Megater  -  Cm3 . 16 

11  Energy  Density  (Energy/ Volume  x  4*R2,  in  Units  of  Megater  - 
Cm2).  Wall  at  8.0  Cm,  Ganna  ■  2.682  Area  Under  Curve 

Equals  Energy  in  Units  of  Megater  -  Cm3  .  17 

TABLES 

Table  Title 

1  Wall  at  6.0  CentlMters;  7  *  2,682  .  18 

2  Wall  at  8.0  Centimeters;  Variable  7 .  20 

3  Wall  at  l6.0  CentlMters;  7  ■  2.682  . 21 

4  Wall  at  24.0  CentlMters;  7  ■  2.682  .  22 

5  Wall  at  32.849  CentlMters;  7  -  2.682  .  23 


ill 


NOLTR  62-19 


INTRODUCTION 

1.  In  carrying  out  a  program  of  investigating  the  characteristics  of 
blast  in  a  rarefied  atmosphere ,  it  is  useful  to  have  available  a  description 
of  the  limiting  case  of  zero  density  in  the  medium  surrounding  the  charge; 
that  is,  an  explosion  in  vacuum.  In  this  case,  the  explosion  consists  only 
of  the  rapidly  expanding  detonation  products,  and  the  loading  produced  upon 
a  target  is  the  result  only  of  the  impacting  detonation  gases,  without  the 
usual  atmospheric  blast  wave.  Numerical  calculations  have  been  carried  out 
at  NOL  on  the  IBM- 704  digital  computer  to  determine  the  characteristics  of 
chemical  explosions  in  vacuum,  and  the  resulting  loadings  on  a  rigid  target. 

2.  The  computational  code  used  is  capable  of  solving  time -dependent 
problems  in  one  space  dimension  with  plane,  cylindrical,  or  spherical 
symmetry.  To  effect  a  compromise  between  the  complex  configurations  occur¬ 
ring  in  cases  of  practical  interest  and  the  capabilities  of  the  machine  code, 
the  following  geometry  was  decided  upon:  A  spherical  charge  of  high  explo¬ 
sive  (in  this  case,  pentolite)  is  detonated  at  its  center,  and  the  detona¬ 
tion  products  expand  outward  with  spherical  symmetry.  At  a  given  distance 
from  the  origin  of  coordinates  (which  is  at  the  center  of  the  charge)  is  a 
rigid,  spherical  wall,  concentric  with  the  surface  of  the  spherical  charge. 
The  quantities  computed  are  the  pressure  and  energy  of  the  gas  as  functions 
of  distance  and  time;  and  the  pressure-time  history  of  the  gas  at  the  rigid 
target. 


INITIAL  CONDITIONS 

3.  The  computation  is  begun  at  the  instant  the  detonation  shock 
reaches  the  surface  of  the  explosive .  It  is  thus  necessary  to  know  the 
conditions  existing  in  the  explosion  products  at  the  completion  of  the 
detonation  process,  since  these  quantities  most  be  used  as  the  Initial 
conditions  for  the  calculation  of  the  expansion.  The  method  used  to  obtain 
these  initial  conditions  is  that  of  0.  I.  Thy  lor,  who  solved  the  problem  of 
determining  the  mechanical  and  thermodynamic  parameters  behind  a  detonation 
shock  by  postulating  a  similarity  solution  where  all  quantities  are  given  as 
functions  of  r/t  (ref.  (a)).  It  should  be  noted  that  the  Taylor  solution 
assumes  that  a  spherical  detonation  wave  may  be  considered  to  be  a  stable 
shock  of  zero  width,  proceeding  radially  outward  with  constant  detonation 
velocity,  and  initiating  chemical  reaction  in  the  material  over  which  it 
sweeps.  Conditions  change  discontinuous ly  across  the  shock,  and  in  particu¬ 
lar,  the  substance  in  front  of  the  shock  is  completely  unreected.  Whereas 
material  just  behind  the  shock  is  completely  reacted.  Refinements  in  this 
detonation  concept,  such  as  considerations  relating  to  the  existence  of  a 
finite  reaction  zone  behind  the  detonation  shock,  are  not  taken  into  account 
in  this  treatment. 
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4.  A  machine  code  that  uses  the  Ideas  of  the  Taylor  treatment  has  been 
written  at  NOL  (ref.  (b)) ,  and  it  is  this  code  that  is  used  to  furnish  the 
initial  conditions  for  the  calculation  of  the  expansion  of  the  explosion 
products . 


EQUATION  OF  STATE 

5.  Many  attests  have  been  made  to  develop  a  satisfactory  equation  of 
state  for  the  burnt  gases  behind  the  detonation  shock,  ranging  from  simple 
perfect  gas  equations  to  elaborate  solid  state  equations.  It  was  felt  that 
excessively  complex  equations  of  state  would  not  be  suitable  for  the  present 
study,  and  Instead,  reliance  has  been  placed  on  certain  results  of  Deal 
(ref.  (c),  (d))  and  considerations  of  Flckett  and  Wood  (ref.  (e))  that 
furnish  some  evidence  to  support  the  belief  that,  for  certain  explosives,  a 
constant-gamma,  ideal  gas  lav  may  be  used  for  detonation  gas  pressures  above 
500  bars.  The  gamma  for  these  gases  is  usually  equal  to  approximately  3. 
Since,  however,  we  are  concerned  here  with  an  expansion  of  the  burnt  gases 
into  a  vacuum,  the  pressure  will  certainly  fall  well  below  500  bars,  and  the 
equation  of  state  must  be  extended  into  this  low  pressure  region. 

6.  One  possibility  that  immediately  suggests  itself  is  simply  to  let 
7  be  constant  over  the  whole  range  of  pressures  encountered  in  the  problem. 
Another  procedure  is  to  take  advantage  of  the  fact  that  at  extremely  low 
pressures,  gases  tend  to  behave  ideally,  with  a  constant  value  of  7  deter¬ 
mined  only  by  the  molecular  structure  of  the  substance.  Since  the  detona¬ 
tion  gases  will  contain,  in  general,  mixtures  of  diatomic  and  triatomlc 
gases,  the  value  of  7  will  be  an  average,  on  a  molar  basis,  of  various  ideal 
gas  gammas.  A  representative  value  for  the  7  of  detonation  gases  at  low 
pressure  is  7  *  1.35*  The  procedure  is,  then,  to  construct  a  function  that 
gives  7  as  a  function  of  pressure,  p,  and  has  the  following  properties  (the 
particular  constants  chosen  refer  to  the  explosive  pentolite): 


7  *  2.682  (for  pentolite)  *  7g  P^Pg  *  1000  bars 


7  “  1.35  »  7X 


p  P^  «  100  bars 


0 


p  ■  1000  bars 


dp 


0 


p  *  100  bars 


r 


2 


4 


71  -  1-35 


72  «  2.682 


100  bars 


1000  bars 


Pig.  1  General  Form  of  7  Versus  p  Curve 
The  section  (A)  of  this  curve  is  fitted  by  a  cubic  equation: 

7  ■  Aq  +  A^p  +  Agp2  ♦  A^p3 


where  the  constants  A^  are  given  by 


7 2  Pj  (Pj  -  3P2)  -  7j  P2  (P2  -  3?!) 
(Pi  -  p2)3 

6PjP2  (y2  -  7X) 

(Pi  -  P2)3 
i(7i  -  72)  (Px  +  P2) 

<ri  -  p2>3 
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2(r2  -  7i) 

(pi  “  p2^3 


7.  The  functional  form  (l)  for  7  -  r(p)  it  thtn  uttd  in  tht  ideal  gat 
equation  of  state  £  M  pT 

where  E  -  Energy/ pi 
p  ■  Pressure 
▼  -  Specific  volume. 

Thus,  two  s cheats  hate  been  suggested  for  construction  of  the  equation  of 
state,  one  which  utilises  a  value  of  7  which  is  constant  over  the  complete 
pressure  range,  and  another  which  fits  the  interval  between  a  constant  high- 
pressure  7  and  a  constant  low-pressure  7  by  naans  of  a  cubic  curve  giving  7 
in  terns  of  pressure.  Most  calculations  have  been  carried  out  using  both 
forns,  so  that  an  indication  nay  be  obtained  of  how  sensitive  various  param¬ 
eters  of  Interest  are  to  the  choice  of  equation  of  state. 


maaraoi  of  as  hmekecal  oaeurmoi 

0.  The  couqputations  were  carried  out  on  an  m-70h  sig<+-mT  computer, 
using  a  revised  version  of  the  "KD-OOIB,"  a  hydrodynamic  code  originally 
developed  at  the  Lawrence  Radiation  Laboratory,  Livermore,  California 
(ref.  (f)).  The  prototype  difference  equations  are  those  of  von  luunaim 
and  Riehtayer,  and  the  technique  used  to  allow  cowutatlon  beyond  shocks  is 
the  von  Veunann  artificial  viscosity  method  (ref.  (g),  (h)).  In  this  msthod, 
a  quantity  analogous  to  a  viscosity  is  introduced  into  the  equations  for  the 
purpose  of  smoothing  discontinuous  shocks  so  that  they  extend  over  a  fixed 
nvmiber  of  spatial  zones  while  meintm-i»i<Mg  the  validity  of  the  Bankine- 
Bugoniot  equations,  which  relate  parameters  on  both  sides  of  the  shock. 
Consequently,  the  machine  program  does  not  find  it  necessary  to  deal  with 
discontinuous  changes  in  the  calculated  quantities,  and  any  shocks  that 
appear  in  the  problem  are  properly  propagated,  although  with  some  rounding 
of  the  normally  steep  shock  front. 

9.  The  calculation  is  carried  out  by  dividing  the  spherical  region 
into  a  given  number  of  spatial  tones,  not  necessarily  all  of  the  same  radial 
length.  At  each  interface  a  mass  particle  is  localised  and  a  given  particle 
will  always  remain  at  the  same  interface.  Each  mass  particle  has  associated 
with  it  a  velocity,  so  that  velocities  are  calculated  at  the  interfaces; 
quantities  like  pressure,  specific  voluns,  and  total  energy  density  are 
localised  in  the  regions  between  interfaces,  that  is,  at  the  centers  of 
spatial  zones. 
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10.  When  the  expanding  gases  reach  the  spherical  wall,  the  Telocity  of 
the  boundary  is  set  equal  to  zero,  and  the  kinetic  energy  thus  lost  io  added 
to  the  last  spatial  zone  as  internal  energy,  thus  raising  the  pressure  in 
this  last  zone.  Hie  pressure  continues  to  rise  as  the  gases  move  toward  the 
wall,  and  a  compression  shock  moves  inward  toward  the  origin. 

11.  The  HOL  computing  facility  includes  equipment  for  the  display  on 
a  cathode  ray  tube  screen  of  various  quantities  as  they  are  computed. 
Photographs  of  such  plots  are  presented  in  this  report,  and  include  distri¬ 
butions  of  pressure  and  energy  density  with  distance  at  various  tines. 


RESULTS 

12.  Before  metrical  data  are  presented,  a  brief  qualitative  descrip¬ 

tion  of  the  results  will  be  given.  Upon  the  completion  of  the  detonation 
process,  the  explosion  gas  expands  into  the  vacua*,  with  a  consequent  drop 
in  pressure  at  the  gas  boundary  and  an  increase  in  the  velocity  of  the 
boundary.  The  pressure  at  the  gas -vacuum  interface  quickly  drops  to  an 
essentially  zero  value,  and  the  boundary  acquires  a  constant  radial  velocity 
As  the  gas  expands,  the  kinetic  energy  in  the  system  increases,  at  the 

expense  of  the  internal  energy.  When  the  gas  boundary  strikes  the  wall,  a 

compression  shock  is  formed  that  proceeds  inward  toward  the  origin.  At  any 
given  time,  the  pressure  between  the  shoek  and  the  wall  appears  to  be  fairly 
constant,  though  this  constant  level  decreases  with  time.  As  the  shock  pro¬ 
ceeds  inward,  most  of  tbs  kinetic  energy  of  the  gas  is  converted  back  into 
internal  energy.  When  the  reflected  shock  reaches  the  origin,  It  is 
reflected  hack  into  the  gas  as  another  compression  shock.  These  shocks  are 
reflected  hack  and  forth  between  the  wall  and  the  origin  until  some  sort  of 

stable  equilibrium  is  reached.  Since  the  entire  process  takes  place  within 

a  fixed,  constant  volume,  it  is  expected  that  the  final  equilibrium  state 

Is  one  of  stationary  gas,  with  energy,  pressure,  and  volume  related  by  the 
ideal  gas  law  having  the  appropriate  value  of  7.  However,  the  computations 
reported  upon  here  do  not  go  far  enough  in  time  to  demonstrate  this  final 
equilibrium. 

13.  The  calculations  made  refer  to  a  one-pound  sphere  of  pentolitc, 
density  >1.65  gm/ cc ,  detonated  at  its  center,  for  various  values  of  wall 
distance.  Pressure  and  energy  density  curves  are  given  in  Figure'  8-11, 
plotted  against  radial  distance  from  the  origin. 

14.  Bach  point  represents  the  value  of  the  quantity  being  considered 
at  a  point  half-way  between  the  positions  of  adjacent  mass  points,  and 
always  remains  associated  with  the  same  two  mass  points.  Thus,  the  plots 
present  not  only  a  pressure  (or  energy)  versus  distance  picture,  but  also 
indicate  the  motion  of  the  mass  particles  in  time.  The .energy  quantity 
plotted  is  an  energy  per  unit  volume,  multiplied  by  k*r  ,  where  r  is  th* 
radial  distance  from  the  origin.  Thus,  the  area  under  the  curve  between 
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any  two  radial  distance  values  is  equal  to  the  energy  contained  in  the  cor¬ 
responding  portion  of  the  system.  The  interchange  between  internal  and 
kinetic  energy  thus  becomes  clearly  visible  from  the  plots. 

15.  One  can  see  from  the  pressure  plots  the  rounding  of  the  shock 
front  and  the  pressure  oscillations  immediately  behind  the  shock  front; 
these  are  characteristics  peculiar  to  the  particular  computational  method 
used,  namely  the  von  Neumann -Richtnyer  ”q"  method. 

16.  In  the  tables,  numerical  results  are  given  for  various  wall 
distances  from  the  origin.  The  quantities  listed  are  time,  in  microseconds, 
with  the  zero  of  the  time  scale  being  taken  as  the  instant  in  which  the 
detonation  shock  reaches  the  outer  surface  of  the  explosive  charge;  the 
kinetic  energy,  internal  energy,  and  total  energy  (in  megabars -cm3)  present 
in  the  system  at  a  given  time;  and  the  shock  position  (in  centimeters)  and 
wall  pressure  (in  bars),  as  functions  of  time.  Due  to  the  finite  difference 
method  of  calculation,  the  total  energy  is  not  exactly  conserved,  as  it 
should  be  in  a  precise  computation.  The  range  of  variation  of  the  total 
energy  can  be  seen  to  be  small  relative  to  the  magnitude  of  the  energy 
transferred  between  the  kinetic  and  internal  modes.  Table  1  and  Table  2 
both  refer  to  a  wall  at  8.0  centimeters,  but  the  calculations  given  in 
Table  2  use  the  alternative  variable  7  equation  of  state.  Tables  4,  and  5 
refer  to  walls  at  the  distances  l6.0  cms.,  24.0  cms.,  and  32.85  cms.,  respec¬ 
tively,  and  each  represents  a  tabulation  using  the  constant  gamma  equation  of 
state.  Figures  2-6  present  plots  of  the  wall  pressure  versus  tine  for  various 
wall  distances,  using  both  forms  of  the  equation  of  state.  For  small  wall 
distances,  where  the  pressure  at  the  wall  remains  well  above  1000  bars,  the 
two  calculations  give  fairly  identical  results.  However,  for  intermediate 
distances,  when  the  wall  pressure  lies  in  the  100-1000  bar  range,  the  results 
begin  to  diverge,  and  for  large  wall  distances,  with  pressures  well  below  100 
bars,  the  variable  7  equation  of  state  gives  substantially  lower  pressures. 
This  behavior  can  be  seen  most  clearly  from  Figure  7,  which  gives  the  peak 
wall  pressure  versus  distance,  for  both  equations  of  state.  The  true  equation 
of  state  probably  yields  curves  which  lie  between  these  two  extremes. 


SCALING 

17.  The  only  characteristic  length  that  appears  in  this  problem  is  the 
initial  charge  radius,  and  therefore  all  lengths  may  be  expressed  in  terms 
of  this  quantity.  Furthermore,  the  initial  distribution  of  mechanical  and 
thermodynamic  parameters  (that  is,  the  Taylor  Wave)  depends  only  on  the 
properties  of  the  solid  explosive  and  is  self -similar  with  respect  to  the 
charge  radius.  Consequently,  the  results  given  in  this  report  for  specific 
distances,  times,  and  explosive  weights  also  apply  to  situations  in  which 
times  and  distances  are  multiplied  by  a  factor  k,  and  the  explosive  weight 
is  multiplied  by  k’.  For  example,  it  can  be  seen  from  Figure  7  that  for  one 
pound  of  pentolite,  a  peak  pressure  of  10  bars  (145  psi)  is  obtained  at  the 
wall  when  the  wall  is  approximately  100  cms.  from  the  origin.  Thus,  a  512- 
pound  pentolite  charge  (k*8)  will  generate  a  pressure  of  145  psi  at  a  wall 
located  approximately  800  cms.  or  26  feet,  from  the  center  of  the  charge. 
Such  a  pressure  is  believed  to  be  high  enough  to  desroy  a  spacecraft. 
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CONCLUDING  REMARKS 

18.  The  results  presented  here  constitute  a  preliminary  Investigation 
of  the  loading  of  structures  by  explosions  in  rarefied  atmospheres.  The 
next  step  is  to  run  calculations  with  air  of  varying  densities  in  the  space 
between  the  explosive  and  vail,  and  to  study  the  air  or  blast  effect.  It 
would  also  be  extremely  interesting  to  substitute  for  an  idealized  rigid 
target,  a  substance  that  may  be  compressed,  or  that  may  be  pushed  as  a 
result  of  the  explosion.  Another  goal  for  the  future  is  the  use  of  more 
practical  geometries,  perhaps  with  the  use  of  2-dimensional  machine  codes, 
and  more  sophisticated  equations  of  state  for  the  detonation  gases,  such  as 
that  of  Klstiakowsky  and  Wilson. 
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FIG.  7  PEAK  PRESSURE  VS  WALL  DISTANCE 
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FIG  9  ENERGY  DENSITY  (ENERGY  /  VOLUME  X  4*R2  IN  UNITS  OF  MEGABAR  - 
CM2).  WALL  AT  8.0  CM,  GAMMA*  2.682. AREA  UNDER  CURVE  EQUALS 
ENERGY  IN  UNITS  OF  MEGABAR -CM?  I-LB  PENTOLITE 
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FIG.  10  ENERGY  DENSITY  (ENERGY/ VOLUME  X  4rR2,  IN  UNITS  OF  MEGABAR  - 
CM2).  WALL  AT  8.0  CM,  GAMMA  -  2.682.  AREA  UNDER  CURVE  EQUALS 
ENERGY  IN  UNITS  OF  MEGABAR- CM3.  I -LB  PENTOLITE 
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Table  1 

WALL  AT  8.0  CENTIMETERS;  7  “  2.682 
ONE  POUND  PEROLITB 


Tine 

(tisecs) 

Kinetic  Internal 

Energy  Energy 

< - (Megabars -Cm3)- 

Total 
Energy 
- > 

Shock 

Position 

(Cms) 

Wall 

Pressure 

(Bars) 

If.  939  (Wall  Struck) 

5-019 

15.745 

5.231 

20.976 

7.559 

3753.9 

6.840 

11.587 

9.602 

21.189 

7.037 

I8a2.0 

7-954 

9.178 

12.007 

21.185 

6.890 

a628.1 

8.951 

7.438 

13.752 

21.190 

6.662 

22603.8 

10.028 

6.004 

15.181 

21.185 

6.400 

22712.4 

11.008 

5.033 

16.147 

a.  180 

6.087 

22096.0 

12.031 

4.242 

16.934 

a.  176 

5.961 

a4oi.o 

13.016 

3-578 

17-597 

21.175 

5.630 

20648.7 

lU.016 

2.925 

18.249 

a.  174 

5.325 

19934.7 

15.031 

2.306 

18.866 

a.  172 

5.220 

19310.1 

17.019 

1.338 

19.829 

a.  167 

4.804 

19062.1 

18.542 

0.825 

20.338 

a.  163 

4.332 

19300.8 

20.022 

0.487 

20.674 

a.  161 

3.994 

19348.7 

21.037 

0.339 

20.819 

a.  158 

3.768 

19074.0 

22.010 

0.277 

20.879 

a.  156 

3-512 

18736.3 

22.518 

0.271 

20.884 

21.155 

3.511 

18642.8 

23.025 

0.260 

20.874 

a.  154 

3.261 

l84a.5 

24.041 

0.330 

20.823 

21.153 

3.002 

18092.6 

25.018 

0.407 

20.745 

a.  152 

2.723 

17797-6 

27.042 

0.619 

20.532 

a.  151 

2.031 

16920.9 

28.021 

0.715 

20.436 

a.151 

I.616 

16660.4 

29.268 

0.787 

20.366 

a.  153 

0.733 

16336.5 

30.250 

0.732 

20.422 

a.  154 

0.191 

16026.6 

First  Shock  Reaches  Center 

31.015 

0.589 

20.568 

21.157 

1.374 

15582.0 

32.286 

0.449 

20.709 

a.  158 

1.885 

15066.6 

33-264 

0.323 

20.835 

a. 158 

2.565 

14469.3 

34.271 

O.232 

20.925 

21.157 

2.826 

14264.4 

35.196 

0.178 

20.979 

21.157 

3.270 

14126.7 

39.169 

0.226 

20.927 

21.153 

4.869 

13308.3 

41.687 

0.461 

20.690 

a.151 

6.018 

13272.5 

44.213 

0.662 

20.491 

21.153 

7-757 

15711.1 

Reflected  Shock  Reaches  Wall 
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Table  1  (Cont'd) 


Kinetic  Internal  Total  Shock  Wall 

Time  Energy  Energy  .  Energy  Position  Pressure 
(usees)  < - ( Megabars  -Cur) - >  (Cms)  (Bars) 


46.701 

0.400 

20.756 

21.156 

6.479 

19734.3 

49.192 

O.183 

20.971 

21.154 

5-571 

19701.0 

51.683 

0.119 

21.035 

21.154 

4.8l6 

18890.0 

54.174 

0.148 

21.005 

21.153 

3-977 

17887.8 

56.712 

O.237 

20.915 

21.152 

3-247 

17052.1 

59.203 

O.368 

20.785 

21.153 

2.179 

16369.2 

Reflected  Shock  Reaches 

Center 

61.693 

0.4l6 

20.738 

21.154 

1.078 

15830.9 

64.171 

0.178 

20.978 

2L.156 

1.674 

15266.4 

66.715 

0.093 

21.061 

21.154 

2.916 

14912.2 

69.206 

0.131 

21.022 

21.153 

4.110 

14606.4 

71.698 

0.217 

20.936 

21.153 

4.995 

14387.0 

74.190 

0.314 

20.840 

21.154 

6.099 

14488.3 

76.681 

0.358 

20.798 

2L.156 

7.766 

17132.0 
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Table  2 

WALL  AT  6.0  CENTIMETERS;  Variable  y 
OKI  FOUXD  PHTPOLITE 


Kinetic 

Internal 

Total 

Shock 

Wall 

Time 

Energy 

Energy  „ 

Energy 

Position 

Pressure 

(tisecs) 

4 - (Megabars -Cm13)- 

- > 

(Cms) 

(Bars) 

4.984  (Wall  Struck) 

5.514 

14.747 

6.499 

21.246 

7.466 

11333.3 

6.507 

12.473 

8.775 

21.249 

7.152 

16231.1 

7-993 

9.187 

12.053 

21.240 

6.902 

21629.6 

8.989 

7.470 

13.776 

21.246 

6.669 

22385.9 

10.015 

6.085 

15.156 

21.241 

6.406 

22476.4 

11.006 

5.085 

16.152 

21.237 

6.096 

22131.2 

15.026 

2.336 

18.892 

21.226 

5.230 

19301.1 

16.014 

1.805 

19.420 

21.225 

4.933 

19006.3 

16.443 

1.603 

19.621 

21.224 

4.877 

19168.2 

19-021 

0.705 

20.513 

21.218 

4.356 

19431.3 

22.026 

0.279 

20.932 

21.211 

3.512 

18704.2 

25.036 

0.4l4 

20.792 

2L.206 

2.716 

17746.6 

25.380 

0.448 

20.759 

21.207 

2.696 

17521.3 
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Table  3 

WALL  AT  16.0  CENTIMETERS  j  7  -  2.682 
ORE  POUND  PENTOUTE 


Tine 

(usees) 

Kinetic  Internal 

Energy  Energy  - 

4 - (Megabars-Cnr) 

Total 

Energy 

l - > 

Shock 

Position 

(Cms) 

Wall 

Pressure 

(Bars) 

14.892  (Expanding  Gases  Strike  Wall) 

16.243 

19.042 

2.277 

a.  319 

14.499 

1075-8 

21.223 

13.817 

7.508 

a. 325 

13.290 

1966.5 

26.262 

9-104 

12.214 

a.  318 

12.397 

2447.3 

33.757 

4.653 

16.653 

a. 306 

11.197 

2506.8 

43.745 

1.091 

20.201 

a.  292 

9.660 

2644.4 

50.727 

0.315 

20.971 

a. 286 

8.231 

2561.4 

55-737 

0.481 

20.800 

a. 281 

7-230 

2395.4 

60.729 

0.950 

20.328 

a. 278 

5-^35 

2a2.7 

65.733 

1.428 

19.852 

a.  280 

3.601 

2056.1 

68.198 

1.546 

19.737 

a. 282 

1.638 

1976.1 

Reflected  Shock  Reaches  Center 

73.305 

1.04i 

20.251 

a. 292 

3.980 

1763.9 

78.209 

0.468 

20.824 

a. 292 

6.i4l 

1663.2 

83.239 

0.210 

21.077 

a.  287 

8.208 

1555.4 

88.270 

0.442 

20.840 

a.  282 

10.085 

1507.0 

93.300 

0.959 

20.323 

a.  282 

12.4912 

1485-7 

98.214 

1.270 

20.022 

a.  292 

15.^53 

2259.3 

Reflected  Shock  Reaches  Wall 

100.785 

I.063 

20.229 

a.  292 

13.138 

2434.2 

105.763 

0.645 

20.646 

a. 291 

11.432 

2507.5 

110.755 

0.323 

20.967 

a.  290 

10.279 

2450.2 

115.747 

0.155 

21.133 

a. 288 

9.087 

2399.5 

120.739 

0.200 

a. 087 

a. 286 

7.760 

2336.4 

125.731 

0.445 

20.840 

a. 285 

5.752 

2231.6 
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Table  4 

HALL  AT  24.0  CENTIMETERS;  7  -  2.662 
OH  POUND  PEMTOLITE 


Time 

(usees) 

Kinetic  Internal 
Energy  Energy  . 

< - (Megabar  8 -Ob'5)- 

Total 
Energy 
- > 

Shock 

Position 

(Cbs) 

Wall 

Pressure 

(Bars) 

24.848 

(Expanding 

;  Oases  Strike 

Vail) 

27.859 

18.861 

2.455 

21.316 

21.785 

265.4 

29.085 

18.114 

3.241 

21.355 

21.201 

379.4 

34.150 

14.668 

6.685 

21.353 

20.158 

515.6 

39.122 

11.468 

9.880 

21.348 

19.005 

636.7 

41.630 

10.051 

11.297 

21.348 

18.748 

673.5 

46.601 

7-687 

13.655 

21.342 

18.081 

714.1 

51.598 

5.819 

15.516 

21.335 

16.734 

726.6 

61.372 

2.760 

18.581 

21.341 

15.601 

731-0 

66.430 

1.630 

19.706 

21.336 

14.826 

773*^ 

71.378 

0.867 

20.465 

21.332 

14.334 

805.6 

76.381 

0.428 

20.900 

21.326 

13.194 

794.9 

81.383 

0.301 

21.022 

21.323 

12.643 

765.7 

91.377 

0.775 

20.541 

21.316 

9.560 

693.9 

101.487 

1.567 

19.749 

21.316 

5.835 

618.4 

106.433 

1.810 

19-511 

21.321 

2.623 

579.6 

108.864 

1.800 

19.522 

21.322 

0.777 

561.2 

Reflected  Shock  Reaches  Center 

111.403 

1.584 

19.743 

21.327 

2.713 

535-2 

113.922 

1.330 

20.003 

21.333 

5.335 

517.8 

115-741 

1.142 

20.192 

21.334 

6.740 

501.3 
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Table  5 

WALL  AT  32.849  CWHT1METBBS;  y  -  2.682 
OHE  POU1TD  PSRTOLCTB 


Kinetic 

Internal 

Total 

Shock 

Wall 

Tine 

Energy 

Energy  . 

Energy 

Position 

Pressure 

< - 

-(Megabars  -Car )— 

- > 

(Cns) 

(Bars) 

35.795 

(Expanding 

Gases  Strike  Wall} 

41.547 

18.153 

3.297 

21.450 

26.662 

139.076 

46.509 

15.823 

5.628 

a.  451 

27.710 

174.486 

51.551 

13.453 

7.995 

a. 448 

26.824 

a5-729 

56.484 

11.291 

10.152 

a. 443 

25.801 

245.894 

61.455 

9-377 

12.065 

a. 442 

24.655 

265.639 

66.503 

7-728 

13.707 

a.  435 

24.461 

275.381 

71.551 

6.339 

15.090 

a.  429 

23.266 

279-206 

76.599 

5.085 

16.342 

a. 42? 

22.362 

280.554 

86.557 

2.822 

18-597 

a.  419 

21.257 

285.344 

91.567 

1.946 

19.470 

a.4i6 

20.673 

298.743 

96.590 

1.257 

20.156 

a.  413 

20.002 

314.907 

101.524 

O.763 

20.647 

a.  4io 

19.197 

318.806 

106.580 

0.440 

20.966 

a.4o6 

18.165 

316. l4l 

111.574 

0.297 

21.106 

a.  403 

17.466 

308.708 

116.569 

0.326 

21.073 

21.399 

16.549 

299.  a2 

121.563 

0.487 

20.910 

21.397 

15.537 

288.742 

129.107 

0.894 

20.501 

21.395 

13.232 

271.972 

141.446 

1.660 

19.734 

a. 394 

8.151 

243.706 

146.455 

1.890 

19.507 

a.  397 

5.856 

232.512 

148.986 

1-953 

19.446 

a.  399 

3.608 

225.779 

151.403 

1.956 

19.445 

a.4oi 

1.143 

220.a3 
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